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ABSTRACT 



We present a multicomponent model to explain the features of the pulsed emission and spectrum of the Crab Pulsar, on the basis of X and y-ray 
observations obtained with BeppoSAX, INTEGRAL and CGRO. This model explains the evolution of the pulse shape and of the phase-resolved 
spectra, ranging from the optical/UV to the GeV energy band, on the assumption that the observed emission is due to more components. The 
first component, Co, is assumed to have the pulsed double-peaked profile observed at the optical frequencies, while the second component, Cx, 
is dominant in the interpeak and second peak phase regions. The spectra of these components are modelled with log-parabolic laws and their 
spectral energy distributions have peak energies at 12.2 and 178 keV, respectively. To explain the properties of the pulsed emission in the MeV- 
GeV band, we introduce two more components, Co y and Cx y , with phase distributions similar to those of Co and Cx and log-parabolic spectra 
with the same curvature but peak energies at about 300 MeV and 2 GeV. This multicomponent model is able to reproduce both the broadband 
phase-resolved spectral behaviour and the changes of the pulse shape with energy. We also propose some possible physical interpretations in 
which Co and Cx are emitted by secondary pairs via a synchrotron mechanism while Coy and Cx y can originate either from Compton scattered 
or primary curvature photons. 

Key words. Stars: neutron - pulsars: individual: Crab pulsar (PSR B0531+21) - X-rays: stars - Gamma-rays: observations 



1. Introduction 

The origin of the high energy emission of rotation-powered 
pulsars is still an unsolved problem. One of the main difficul- 
ties is related to the description of the phase and energy distri- 
butions of the pulsed signal, which depends on both a physical 
and geometrical modelling of the magnetosphere. It is still un- 
clear, for example, whether electrons (or positrons) are acceler- 
ated and radiate streaming out the polar cap regions (Ruderman 
& Sutherland 1975, Salvati & Massaro 1978, Sturner & Dermer 
1994, Daugherty & Harding 1994, 1996; Muslimov & Harding 
2003) or in the outer gaps (Cheng, Ho & Ruderman 1986a,b; 
Chiang & Romani 1994, Romani & Yadigaroglu 1995, Cheng, 
Ruderman & Zhang 2000; Zhang & Cheng 2002, Hirotani, 
Harding & Shibata 2003). Another important problem con- 
cerns how quantum processes, like magnetic pair production 
(Erber 1966) and photon splitting JAdler 197U . modify the 
high energy y-ray spectrum. 

The Crab pulsar (PSR B053 1 +2 1) has been the best studied 
object of this class since its discovery (Staelin & Reifenstein 
1968) and the amount of data collected is rich enough to search 
for a detailed physical picture of its emission properties. It is 
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well known that the pulse shape of Crab has a characteristic 
double peak structure, with a phase separation of 0.4, detected 
from the radio band to y rays and changing with energy. A very 
remarkable feature is that the so-called second peak, hereafter 
P2, in the X and soft y-ray ranges becomes progressively higher 
than the first peak (PI). A similar increase is also evident in 
the emission between the peaks, usually named the Interpeak 
region (Ip) or bridge (see Fig. 1). Above about 10 MeV, PI is 
again the dominant feature. A satisfactory explanation for these 
changes has not been found so far. 



On the basis of high quality BeppoSAX data, covering a 
wide energy range from 0. 1 to about 300 keV, we proposed a 
two component model (Massaro et al. 2000, hereafter MCLM) 
to interpret this behaviour. In the same paper we studied in de- 
tail the energy spectrum of the core of PI (corresponding to a 
phase interval having a width of only 0.027 around the maxi- 
mum) which shows a continuous steepening at high energies. 
We found that in the energy range 0.1-300 keV this spectral 
distribution is well represented by a parabolic law in a dou- 
ble logarithmic plot (hereafter log-parabola) with a rather mild 
curvature. The extrapolation of this model in the y-ray range, 
however, fails to reproduce the data and a more complex mod- 
elling is required. 
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Fig. 1. The pulse profile and the phase dependent photon index of the Crab pulsar observed with the four Narrow Field Instruments 
of BeppoSAX and with ISGRI-INTEGRAL experiment. The respective energy ranges are indicated in the upper panels. All the 
profiles are normalized to unity at the maximum of PI. Note the change of relative intensity of P2 and Ip with respect to PI and 
the increase of the photon indices. 



In a subsequent paper Kuiper et al. (2001) introduced three 
components to describe the spectrum up to the y-ray data ob- 
tained with COMPTEL and EGRET on board ComptonGRO. 
These authors based the analysis mainly on finding best fits of 
the spectral distributions in rather narrow phase intervals and 
found that the spectral variation of the pulsed emission with 
phase can be modelled by two log-parabolic components with 
a relative normalization changing with phase. A further third 
power law component, having a photon index equal to 2.07, 
was necessary for the emission of PI and P2 to to match the 
EGRET data at energies higher than about 10 MeV. 
In this paper we develop a model able to describe the phase and 
spectral distributions of the emission over a frequency interval 
from the optical frequencies to GeV range. This model is an ex- 
tension of that presented in MCLM and it based on the results 
of a new detailed data analysis of many BeppoSAX observa- 
tions which includes PDS data from March 1999 to April 2001 
not considered by MCLM. The timing accuracy has been veri- 
fied using detailed pulse profiles obtained from RXTE archive 
data. Moreover, to extend the energy range, we analysed several 
more recent observations performed with the IBIS-ISGRI ex- 
periment on board the INTEGRAL satellite and considered the 
results of Kuiper et al. (2001) on the COMPTEL and EGRET 
observations up to a few GeV. The main goal of our work is 



the definiton of a scenario that can be used to develop more 
detailed physical models of the Crab pulsar high-energy emis- 
sion. 



2. Observations and data reduction 

2.1. BeppoSAX 

BeppoSAX observed the Crab Nebula and its pulsar on several 
occasions during its lifetime, because this very bright source 
was used to monitor the responses of the four Narrow Field 
Instruments (NFI). All the data collected up to 1998 were used 
in MCLM and a complete description of their reduction is pre- 
sented in that work. In the new analysis presented in this pa- 
per we did not consider some early PDS observations, partic- 
ularly those performed on 31 August 1996, during the Science 
Verification Phase, and on 6 April 1998 because the detector 
gain was found to be not precisely set, whereas additional PDS 
data were obtained from other observations performed between 
1999 and 2001, corresponding to an increase of the exposure 
time of about 12% with respect to MCLM. The complete log 
of all BeppoSAX pointings used in the present work is given in 
Table 1. 
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For the imaging instruments we selected all the events 
within circular regions centered on the source position and hav- 
ing radii of 4' (MECS) and 8' (LECS). This choice corresponds 
to using about 90% of the total source signal in both instru- 
ments, but it allows us to apply the best tested spectral response 
matrices. 

Phase histograms of the Crab pulsar were evaluated for 
each NFI and each pointing using the period folding technique. 
The UTC arrival times of all selected events were converted to 
the Solar System Barycentre with the DE200 ephemeris. The 
values of P and P for each observation epoch were derived 
from the Jodrell Bank Crab Pulsar Monthly Ephemeris 1 . We 
constructed a large set of 300 bin phase histograms for each 
energy channel of each NFI. The zero phase was fixed at the 
centre of the first peak, evaluated by Gaussian fits. All the his- 
tograms for the same energy channel of each NFI were then 
added. Before this operation we verified that all the profiles 
of the various observation epochs had fully compatible shapes 
and therefore similar folding accuracy. A summary of all these 
phase histograms is shown in the upper panels of Fig. 1. 

A relevant point we stress when working with spectra ob- 
tained with several instruments is that a proper evaluation of the 
inter-calibration factors is required. For the BeppoSAX NFIs, 
the accurate ground and in-flight calibrations were used to es- 
tablish the admissible ranges for the factors between the MECS 
and the other three instruments: k.ML for the LECS, Icmh for the 
HPGSPC and k MP for the PDS (Fiore et al. 1999). When per- 
forming the spectral analysis, we left these three parameters 
free and found the values £^£=0.74 and ^#=0.88, within the 
ranges given by Fiore et al. (1999). For point sources these au- 
thors gave 0.77 < k M p <0.93, reduced to 0.86±0.03 for sources 
with a PDS count rate higher than 2 ct/s. We found in our best 
fits k M p=0.&2, very close to the above prescription, and that we 
considered fully satisfactory because it was derived with the 
log-parabolic model instead of the simple power law used by 
Fiore etal. (1999). 

2.2. RossiXTE 

To verify the accuracy of BeppoSAX timing, we compared the 
pulse profiles as observed by MECS with that obtained by the 
PCA instrument (Jahoda et al. 1996) on board of the Rossi X- 
Ray Timing Explorer. The latter instrument has a timing ac- 
curacy of 2 yus (Rots et al. 2004) with respect to the space- 
craft clock and an absolute time accuracy of 5-8 /zs with re- 
spect to UTC. We used the RXTE observation of Crab per- 
formed on December 3, 2004. PCA data were obtained in an 
"event mode" that provided 128 energy channels with 250 
jus time resolution. PCA data were screened following stan- 
dard selection criteria excluding time intervals corresponding 
to South Atlantic Anomaly passages, Earth limb closer than 
10 degrees and with angular distance between the source and 
pointing direction larger than 0°.02. The total exposure time 
after all screening criteria was 709 s. After the conversion of 
the arrival times to the Solar System Barycentre by using the 
DE200 ephemeris, phase histograms were obtained by folding 

1 http://www.jb.man.ac.uk/. 



Table 1. Log of BeppoSAX NFIs Crab Pulsar pointing epochs 
and net exposure times. 



Observation 




Exposure Times (s) 




Date 


LECS 


MECS" 


HPGSPC 


PDS 


31 Aug 1996 






28 384 




6 Sep 1996 


5 733 


33 482 






30 Sep 1996 




7 874 




3596 


12 Apr 1997 


1257 


18551 


9 344 


9 338 


8 Oct 1997 


12 096 


30769 


9158 


19 673 


6 Apr 1998 


8 730 


28 694 


12 732 




13 Oct 1998 








17 476 


9 Mar 1999 








13 674 


25 Sep 1999 








15417 


10 Apr 2000 








15 602 


4 Apr 2001 








7 597 


Total exposure 


27 816 


119 370 


59618 


102 373 



" MECS operated with two detectors after May 1997. 



Table 2. Log of INTEGRAL-ISGRI observations of the Crab 
and of net exposure times. 



Start - 


Stop day 


Exposure Times (s) 


07 Feb 2003 


- 09 Feb 2003 


142 300 


10 Feb 2003 


- 11 Feb 2003 


56000 


16 Feb 2003 


- 18 Feb 2003 


61600 


19 Feb 2003 


- 21 Feb 2003 


30700 


22 Feb 2003 


- 24 Feb 2003 


9 800 


25 Feb 2003 


- 26 Feb 2003 


2100 


15 Aug 2003 


- 17 Aug 2003 


37 900 


17 Aug 2003 


- 17 Aug 2003 


21600 



Total exposure 362000 



data with the Jodrell Bank Crab Pulsar Monthly Ephemeris as 
before. 

Fig. 2 shows the RossiXTE-PCA pulse profile, extracted 
in the 2.6-4.1 keV energy band and with 300 phase bins (110 
//s), together with that obtained in the same energy range with 
MECS. Both profiles have been normalized to unity at the max- 
imum of PI, after subtraction of the off -pulse constant level 
(0.6-0.83) and arbitrarily shifted to set the zero phase at the 
centre of the first peak. The comparison shows that the peak 
widths between the two instruments are very similar. The two 
profiles have differences of only a few percent in amplitude 
(Fig. 2, lower panel) and there is no evidence of a systematic 
effect on the timing in the MECS data. We conclude that the 
BeppoSAX profiles accumulated over multiple observations do 
not undergo any significant broadening due either to the clock 
time assignment or phase misalignment among the different 
histograms added together. 

2.3. INTEGRAL 

To gain more information at energies higher than 200 keV 
we considered also the recent data obtained with IBIS-ISGRI 
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Fig. 2. The comparison of the normalized pulse profiles in the 
2.6-4.1 keV energy range of BeppoSAX-MECS (solid his- 
togram) and RossiXTE-PCA (dashed histogram) (upper panel). 
The difference plot in the lower panel has an amplitude of only 
a few percent and no sign of systematic effects on the timing. In 
the upper panel we also show the phase intervals of PI, Ip and 
P2 used in our spectral analysis (solid vertical lines - black) 
and according to Kuiper et al. (2001) (dotted vertical lines - 
red) used to compute the broad-band spectral distributions and 
the P2/P1 and Ip/Pl ratios (see Sect. 4). 

(Lebrun et al. 2003) on board the INTEGRAL satellite. In our 
analysis we considered only the IBIS observations having an 
off-axis angle not larger than 1 which are listed in Table |3 
Phase histograms with 100 bins were obtained for each point- 
ing using, as before, the period folding technique (see Fig. 1). 
The accuracy of INTEGRAL timing with respect to RXTE 
has been verified by Kuiper et al. (2003) The zero phase was 
shifted at the centre of the first peak, coherently with the as- 
sumption for the other instruments. A more complete spectral 
analysis of these data, together with those obtained using other 
INTEGRAL instruments, can be found in Mineo et al. (2006). 

3. The interstellar absorption in the Crab direction 

3.1. X-ray absorption 

An important piece of information for the study of the intrin- 
sic shape of the soft X-ray spectral distribution is the value of 
the equivalent hydrogen column density and the composition of 
the interstellar matter in the Crab direction. In MCLM we as- 
sumed jVh=3.2x10 21 cirT 2 and a mean solar composition. This 
value is higher than those considered in previous papers, usu- 
ally lower than 3xl0 21 cirT 2 (see, for instance, Toor & Seward 
1977, Ride & Walker 1977, Pravdo & Serlemitsos 1981) and 
in agreement with the analysis of ASCA data by Fukuzawa 



PI (First peak) 
Ip (Interpeak) 
P2 (Second peak) 
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Fig. 3. The energy dependence of the photon indices of PI, Ip 
and P2, evaluated in the ranges of the LECS, MECS , HPGSPC, 
two PDS bands (15-80, 80-200 keV) and INTEGRAL-ISGRI 
(50-500 keV, large symbols). The off-pulse photon index, rep- 
resentative of the nebular emission, is also shown. Best-fit lin- 
ear interpolations of BeppoSAX photon indices of PI, Ip and 
P2 are plotted. 



et al. (1997) and with the result of Sollerman et al. (2000) 
who estimated N H = 3.2 ± 0.5 x 10 21 cm 2 from STIS/HST 
observations of the Lyo- line. This column density, however, 
was not confirmed by subsequent studies that indicated an even 
larger A^h and different chemical abundances. Willingale et al. 
(2001) on the basis of EPIC-MOS data found a column den- 
sity of 3.45xl0 21 cirT 2 with an underabundance of oxygen and 
iron of 63% with respect to the solar composition. Weisskopf et 
al. (2004) obtained N H = 4.2xl0 21 crrT 2 with the abundances 
given by Wilms et al. (2000) and using a photoelectric absorp- 
tion model that includes the scattering due to interstellar grains. 

In order to use a more precise estimate of A^h and the 
most suited chemical composition we performed a new analy- 
sis of the X-ray spectrum of the DC emission of Crab (pulse 
phases 0.60-0.83) using LECS and MECS data combined. 
Background, estimated from archive blank fields, was sub- 
tracted from the events. We adopted a simple power law 
(F(E) = K{E I EqY t ) as the best fitting model because in this 
band there is no evidence of a change of the nebular spectral in- 
dex with energy (see next Section). The considered abundances 
and photoelectric cross sections were those of the models 
aneb+wabs (Anders & Ebihara 1982, Morrison & McCammon 
1983) and angr+bcmc (Anders & Grevesse 1989, Balucinska- 
Church & McCammon 1992, vphabs) all implemented in the 
XSPEC package. The results are reported in Table [3] where 
for each model we give the best fit values of A^h and of the 
photon index F, the fraction of oxygen and iron with respect 
to the default value and the reduced^ 2 . All the fits were ac- 
ceptable (reduced x 2 < 1 ) and the column densities were sys- 
tematically higher than 3.2xl0 21 cirT 2 . Photon index values 
were found to be very stable with a mean value of 2.124. We 
also found (see Sect. 6.2) that the residuals at energies lower 



E. Massaro et al.: A multicomponent model of the Crab pulsar 



5 



than 0.8 keV had a systematic deviation which disappeared us- 
ing column density values around 3.6xl0 21 crrT 2 . We decided, 
therefore, to adopt in our analysis this value together with 
the angr+bcmc model, as indicated by other instruments with 
higher spectral capabilities. 

3.2. Optical-UV extinction 

The high estimate of the equivalent column density and the 
possible difference of the chemical composition with respect 
to the solar one could also affect the optical-UV extinction in 
the direction of Crab. The analysis of Sollerman et al. (2000) 
converges around a value of E(B — V) = 0.50 mag and indi- 
cates that the extinction curve closely follows the standard one. 
Thus the N H /E(B - V) ratio is equal to 6.9xl0 21 cirT 2 mag -1 , 
about 30% higher than the value reported by Predehl & Schmitt 
(1995). 

Such a relevant optical-UV intrinsic extinction makes the 
evaluation of the pulsed flux and spectrum quite difficult at 
these frequencies. Several authors (e.g. Percival et al. 1993, 
Sollerman et al. 2000) agree that the dereddened spectrum is 
a power law with a spectral index close to zero. The difference 
to uncorrected data is quite large: in fact, the observed spec- 
tral index is around -2, and small differences in the values of 
E(B - V) or Ay can modify the dereddened slope and flux. Of 
course, the extrapolation of UV data in the X-ray range may be 
affected by large uncertainties. 

4. Phase resolved X-ray spectroscopy 

In MCLM we showed that a simple power law does not pro- 
vide an acceptable representation of the spectral distribution of 
the central bins of PI, corresponding to the narrow phase in- 
terval 0.99-0.01667. This is not a specific property of PI but 
can be extended to the entire pulsed signal as apparent from 
the phase resolved spectroscopy of Fig. [2 where pulse profiles 
and best fit photon indices evaluated in the energy bands of the 
fourNFIs of BeppoSAX and of ISGRI-INTEGRAL are shown. 
Upper panels show the pulse profiles while lower panels show 
the phase evolution of the photon indices: note the well known 
increase of Ip and P2 with respect to PI and that photon indices 
are also increasing with energy: that of PI changes from about 
1.6 to 2.4 and that of the Ip from 1.3 to 1.95. The power law 
fit was applied to the events after the subtraction of a DC level, 
estimated from the phase interval 0.60-0.83, and assumed to be 
representative of the nebular continuum. High resolution X-ray 
images obtained with Chandra (Tennant et al. 2001) showed 
that there is a weak emission from the pulsar also in this phase 
interval. Its mean intensity is of the order of few percent than 
that of the peaks and becomes negligible compared to the neb- 
ular emission when considering data extracted from wider an- 
gular beams, as in our case. 

To study the behaviour of the photon index with energy in 
detail we considered the following three phase intervals rep- 
resentative of PI (0.9833-0.0167), the Interpeak (Ip) region 
(0.14-0.25) and P2 (0.3833-0.4137). This choice was made to 
have intervals wide enough to obtain estimates of F with small 
statistical uncertainties and sufficiently narrow to consider the 
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Fig. 4. The shape of PI in three energy ranges. Maxima have 
been normalized to unity for a better comparison. Note the 
broadening of the peak with a marked asymmetry on the right 
side. Solid lines represent the PI profiles derived using the two 
component model (see Sect. 6.3). 

emission of PI mainly due to one of the components introduced 
by MCLM. In subsequent analysis, when comparing the results 
of our multicomponent model with literature data (see Sect. 
7), it was necessary to consider for PI and P2 phase intervals 
broader than those given above and we adopted the same def- 
inition given by Kuiper et al. (2001). The phase interval of Ip 
was the same. Fig. 2 shows the differences between these two 
sets. 

The resulting values of F in the narrow phase intervals 
are shown in Fig. [3] together with those of the off-pulse, as- 
sumed originating in the nebula. The values of the PI pho- 
ton index vary from 1.72+0.02 in the O.l^l- keV range to 
2.30+0.05 at (80-200 keV); correspondingly, those of Ip vary 
from 1.31+0.10 to 1.98±0. 10 and those of P2 from 1.51+0.03 
to 2.14+0.05. We considered also INTEGRAL-ISGRI data in 
the rather broad interval 50-500 keV and obtained photon in- 
dices fully consistent with the increasing trend of BeppoSAX 
values. The nebular spectrum has a remarkably constant value 
equal to 2.131+0.002, up to the PDS range where it slightly in- 
creases to 2.172+0.006 (15-80 keV) and to 2.23+0.02 (80-200 
keV). Note that the last value agrees very well with that ob- 
tained by Kuiper et al. (2001) from COMPTEL data between 
0.75 and 30 MeV. The very good consistency of the nebular 
photon index measured in LECS, MECS and HPGSPC is an 
indication of the high accuracy of their calibration and of the 
correct choice of the intercalibration constants. 

It is clear from these results that a law more complex than 
a simple power law must be used to describe the broad-band 
spectral properties of the Crab pulsar when the data of all the 
NFIs are considered. A simple power law must be rejcted be- 
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Table 3. Estimates of the equivalent hydrogen column density 7V H , from combined LECS (0. 1—4 keV) and MECS (1 .6-10 keV) 
observations with various absorption models and chemical abundances. 



Model Nun r Fo F Fe X ,- 



aneb+wabs 


3.54 d 


= 0.01 


2.124 = 


t 0.001 










0.89 


angr+bcmc 


3.33 d 


: 0.01 


2.124 = 


t 0.001 










0.99 


angr+bcmc 


3.62 = 


= 0.04 


2.125 = 


t 0.001 


0.81 = 


t0.02 






0.90 


angr+bcmc 


3.63 d 


= 0.04 


2.124 = 


t 0.001 


0.83 = 


t0.03 


0.87 


± 0.05 


0.88 


angr+bcmc 


3.42 = 


= 0.02 


2.123 = 


t 0.001 






0.74 


± 0.04 


0.93 



(*) In units of 10 21 cm" 2 . 



cause of the high values of the reduced^ 2 , found always larger 
than 2. 

In the MCLM analysis of the PI spectrum, good spectral 
fits are obtained using a second order law in the double log 
representation: 



F(E) = A (E/E ) 



-{a+b Log(£/£ )) 



(1) 



where Eq is taken equal to 1 keV, and therefore a corresponds 
to the photon index at this energy, while b measures the curva- 
ture of the spectral distribution. We can also define an energy 
dependent photon index: 



HE) 



dLogF(E) 
dLogE 



a + 2b Log(E/E Q ) 



(2) 



This law implies then a linear relation between F and LogE 
like that apparent for the best fit lines plotted in Fig. [3] which 
have remarkably similar slopes corresponding to a curvature 
parameter b between 0.13 and 0.16, the former corresponding 
to P2 and the latter to Ip. 

We are also interested in the spectral energy distribution, 
defined as E 2 F(E), which for a log-parabolic law has a maxi- 
mum at the energy: 



E p = \Q( 2 - a '» 2b . 



(3) 



Mineo et al. (2006) used a log-parabola for the phase re- 
solved spectroscopy of the INTEGRAL observations of Crab 
and found a significant curvature which was found stable across 
the entire phase interval of the pulsed signal with a mean value 
of 0. 14±0.02. We apply this law to each of the two components 
introduced by MCLM, as described in detail in Sect. 6.2. 

5. Shape and spectrum of the first peak 

Fig. 0] shows the normalized profiles of PI in the three energy 
bands 1.6-3 keV, 16-32 keV and 87-300 keV. As noticed by 
Mineo et al. (1997), there is an evident broadening which in- 
creases with energy. In the lowest energy range the peak pro- 
file is asymmetric with a leading side higher than the trailing 
one, as observed at optical frequencies (see, e.g., Kanbach et al. 
2003). This asymmetry changes in the hard X-rays as the trail- 
ing side becomes progressively higher. The broadening, how- 
ever, can be safely detected also on the leading side. We verified 
that this effect is also present to the same extent in the RXTE 
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Fig. 5. Fine phase resolved photon index across PI in the 1.6- 
10 keV (MECS) and 20-100 keV (PDS) energy ranges. Solid 
lines are the photon indices derived using the two-component 
model as explained in Sect 6.4. 



data and therefore it cannot be due to possible phase misaligne- 
ments of the various BeppoSAX data sets. 

Pravdo, Angelini & Harding (1997), using a RXTE obser- 
vation of the Crab pulsar, found that the central part of PI has 
a spectrum softer than the wings. In MCLM we confirmed 
this result from MECS data, but it was after questioned by 
Vivekanand (2002), who in a new analysis of the RXTE data in 
the range 5-60 keV obtained a stable value of F across PL We 
performed a detailed analysis with a finer phase resolution and 
confirmed the previous results: the photon index in the MECS 
range changes from ~ 1.7 in the wings of PI to 1.82 in the core, 
as shown in the plot of Fig. [5] The same effect is also apparent 
in the PDS data. As it will be clearer in the following section 
this result is relevant to explain the total pulsed emission as due 
to the superposition of two components with different spectra. 
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Phase 

Fig. 6. Upper panel: phase distributions of the components Co 
and Cx with the proper normalizations. Co is derived from 
OPTIMA lightcurve (Kanbach et al. 2003). The four vertical 
bars indicate the phases from fa to fa used in the analytical 
model of Cx- Lower panel: comparison between the observed 
pulse profile in the energy interval (8.0-10.0) keV (MECS 
data) and that obtained using the two component model com- 
puted for the mean energy of 8.85 keV. 



6. The two component model 

To explain at the same time the change of the pulse profile and 
the phase evolution of the photon index, in MCLM we pro- 
posed a simple model in which the observed pulsed emission 
is considered as the superposition of two components having 
different phase distributions and energy spectra. The first com- 
ponent was assumed to have the same pulsed profile observed 
at optical frequencies (Smith et al. 1988), with PI much more 
prominent than P2 and a very low intensity in the Ip region. 
The reason for this choice is that a similar shape is also ob- 
served at energies higher than 30 MeV, despite the profile of 
P2 being slightly different (Fierro et al. 1998). This component 
was called 'optical' (shortly Co)- The other 'X-ray' component 
Cx reaches the greatest relative intensity in the hard X to low- 
energy y-rays and it is necessary to explain the change of the 
P2/P1 ratio. We derived its phase profile to reproduce the ob- 
served one when it is added to Co'- we found that Cx increases 
monotonically up to the phase 0.4 and then it has a rather sharp 
cut-off (Fig. 6). The main properties of the latter component 
were then estimated from the BeppoSAX data, with a fitting 
procedure of several pulse profiles at different energies, as ex- 
plained in detail in MCLM. 

In this new analysis we followed the same approach and 
considered again the same two components with only slight 
changes, as explained in the following subsections. 



6.1. The components' shapes 

We obtained the phase profile of Co from the recent light curve 
observed with the high speed photo-polarimeter OPTIMA with 
a time binning of 112 ps (Kanbach et al. 2003), comparable 
to that used for the BeppoSAX profiles. We used a high res- 
olution digital scanner to convert the OPTIMA data and then 
we applied a rebinning to have the final phase profile with 300 
bins, the same used for the X-ray data. This profile is shown in 
the upper panel of Fig. 6. 

The pulse shape of Cx can be obtained by the difference 
between the observed profiles and Co - However, using the same 
approach as in MCLM we preferred to model Cx by means of 
an analytical expression to make simpler computations. Cx was 
factorized in two functions: 

F x (E,f) = Y(E)g(f,E) , (4) 

where Y(E) is assumed to depend only on the photon energy E 
and represents the spectral distribution of the component to be 
derived from the data as in MCLM, while g(f, E) describes the 
phase behaviour and depends on both variables. In MCLM we 
found that the assumption that g is only a function of the phase 
was not fully satisfied and a moderate change of the phase pro- 
file with the energy will be considered in this paper for a more 
precise description of the data. The whole phase interval of Cx 
is denoted by (fufa) and it is divided by two inner points at the 
phases fa and fa into three segments: in (fa , fa) the Cx shape is 
given by a power law, in (fa, fa) by an exponential function 
joining the first one at fa, and finally a linearly descending 
branch from fa to fa that connects the maximum to the zero 
level of the off-pulse. We have then: 

g(faE) = exp{p(E)(fa-fa)}[j^-^ , (5) 
for fa < f < fa, and 

g(f,E) = exv{p(E)(f-fa)}, (6) 

for/2 </ <fa. 

We emphasize that Cx is non zero also in the phase inter- 
val of PI: this choice is necessary to explain the broadening of 
this peak with increasing energy and the softer spectrum of the 
peak centre with respect to both wings (see Sect. 5.1). Note, 
however, that because of the dominance of PI with respect to 
Ip in Co, is it difficult to establish the actual shape of Cx across 
PI. Our assumption that it has the smooth power law shape of 
Eq. (5) must be considered as one of the simplest approxima- 
tions. 

In MCLM we found that the phase interval boundaries, 
equal to fa = -0.15, fa = 0.21, fa = 0.40, fa = 0.43 and 
the power law exponent s = 0.4 can be considered constant 
with the energy, whereas a better agreement with the data is 
obtained if p is allowed to vary, changing from about 1 1 in the 
keV range to 8 above 150 keV. A simple regression gives the 
relation: 

p(E) = -1.5811 ln(LogE) + 9.318 (7) 
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Fig. 7. Spectral best fit of the four BeppoS AX-NFI data sets of 
Ip with the log-parabolic law of Eq.(9). 



that we used in our calculations. An example of how these two 
components combine to reproduce the observed profile in the 
energy range 8-10 keV is shown in Fig. [6] 

Even though this model has some minor defects of accu- 
racy, it gives an acceptable picture of the broad-band spectral 
and phase properties of the pulsed emission of Crab, as we will 
show in the following sections. This model is not only a useful 
tool to represent a limited data set and it can be successfully 
used to extrapolate the pulse shape and spectra outside the en- 
ergy range in which it is established. Note that Kuiper et al. 
(2001) derived a seven-bin phase profile of their component 
necessary to explain the spectral difference between Ip and the 
peaks and found a structure similar to Cx- 

6.2. Spectral properties of the two components 

To model the phase and spectral properties of the total emission 
as the superposition of Co and Cx we had to estimate their 
spectral distributions. The low intensity of Co in the Ip interval 
can be used to represent its X-ray flux as due only to the Cx, 
while in the first and second peaks both components contribute 
to the total emission, although with different relative intensities. 
Assuming that the spectrum of each component is given by a 
log-parabolic law, we have for PI, Ip and P2: 



F(E) Pl = (A ) P1 E 



-(£i +£ Log£) 



+ (Ax)piE 



-(ax+bx~LogE) 



(8) 



F(E) Ip = (A x ), p E- (aMo ^ (9) 

F(E)p2 = (Ao)p2E- (ao+hohosE) + (Ax)p2E- (ax+hxhosE) (10) 

The parameters to be determined are thus cio, bo, ax and bx, 
besides the normalization factors {Ao)p\, {Ao)pi and (Ax)i p . 
The constants {Ax)p\ and (Ax)p2 can be obtained from (Ax)i p 
and from the shape g(f,E) of Cx- Let Fx(E, f) be the flux of 
Cx component; we have: 



f Ip F x (E,f)df _ (Ax)lp 
J plp2 F x (E,f)df (A x ) Pl ,P2 



(ID 
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Fig. 8. Spectral best fit of the four BeppoS AX-NFI data sets of 
PI with the combined log-parabolic model of Eq.(8) and with 
the UV constraint of Eq.(12). 




channel energy (keV) 

Fig. 9. Spectral best fit of the four BeppoS AX-NFI data sets of 
P2 with the combined log -parabolic model of Eq.(10) and with 
the UV constraint of Eq.(13). 



With the factorization of Eq. (4), we obtain: 



{A x )p\,P2 = (A x )i P 



$ PhP2 g(f,E)df 
I Ip 8(f,E)df 



(12) 



Then, integrating g(f, E) in the proper phase intervals of Ip, 
PI and P2, and obtaining (A x )i p from a spectral fit, we can 
evaluate (A x )p\ and (A x )p2- 

To obtain a fully consistent picture, in addition to the X-ray 
data, the spectral distribution of Co must satisfy the optical-UV 
observations. Sollerman et al. (2000), from STIS/HST observa- 
tions, found that the UV spectrum of the Crab pulsar is remark- 
ably flat for both peaks and can be described by a single power 
law, with an energy index auv = -0. 1 1 . To better constrain our 
spectral parameters we impose the condition that the extrapo- 
lation of the Co log-parabolic spectral distribution in the UV 
range would be consistent with this value and from Eq. (2) we 
derived the following relation between cio and bo'- 



b = 0.21 a - 0.187 



(13) 
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Table 4. Parameters of the two-component, log-parabolic spectral model fits of BeppoSAX data (LECS, 0.1—4.0 keV; MECS 
1.6-10.0 keV; HPGSPC 10-34 keV; PDS 15-300 keV). 



Parameter 


First peak (PI) 


Interpeak (Ip) 


Second peak (P2) 


a 


1.652 ± 0.003 




1.670 ± 0.008 


b 


0.160 ±0.001 




0.164 ±0.001 


a x 




1.28 ± 0.04 




b x 




0.16 ±0.01 




AoO 


0.1660 ±0.0005 




(4.72 ± 0.06) • 10~ 2 


AxO 




(2.3 ±0.1)- 10~ 2 




xl 


1.02 (867 d.o.f) 


0.99 (866 d.o.f) 


1.09 (867 d.o.f) 



(*) photons cm 2 s 1 keV 1 



which was included in the spectral best fitting procedure. The 
resulting best fit spectra with the model of Eqs. (8)-(10) and the 
residuals for Ip, PI and P2 of all NFI data sets of BeppoSAX 
are shown in Figs.0|8j|9] respectively, and the best fit values 
are reported in Table |4] 

The condition of Eq. (13) should not be considered very 
stringent in the evaluation of the spectral parameters. In fact, 
when it is released and only X-ray data are used, the best fit 
values of a and b are consistent with the ones in in Table |4] 
Eq.(13), therefore, must be considered as a useful tool to re- 
duce the parameter space and to avoid possible inconsistencies 
between the optical-UV and X-ray data. Moreover, the actual 
spectral distribution at frequencies much lower than that of the 
peak can deviate from a pure log-parabola and can approximate 
a less curved law like a simple power law. Unfortunately, in this 
frequency range extinction effects are important (see Sect. 3.2) 
and it is not simple to reach a satisfactory picture. 

Three other Crab-like pulsars show curved X-ray spectra 
that can be fitted by log-parabolae with curvature parameters 
very similar to that of Crab: PSR B 1509-58 (b = 0.16±0.04, 
Cusumano et al. 2001), PSR B0540-69 (b=0. 143+0.003, 
de Plaa, Kuiper & Hermsen 2003) and PSR J0537-6910 
(b=0. 15±0.05, Mineo, Cusumano & Massaro 2004). This find- 
ing can be considered an indication that the same radiation 
mechanism is at work in these sources. 

We recall that there are important differences between our 
modelling and the one of Kuiper et al. (2001). These authors 
assumed that the pulsed signal observed above ~30 MeV, mod- 
eled by a single power law, extrapolates down to the X-ray 
range and therefore fitted the phase resolved spectra by means 
of this law plus two log-parabolae. The estimate of the cur- 
vature parameters is thus affected by the power law, and their 
values (0.299 and 0.084) are quite different from those found 
by us and cannot be directly related to the slopes of the linear 
best fits of Fig. 3. The normalization factors of these compo- 
nents are established independently in each phase interval con- 
sidered, increasing the number of free parameters, whereas we 
evaluated the normalization factors of Co and Cx from their 
phase distributions, as explained above. 



6.3. The X-ray shape and the spectral index of P1 

The two componenent model is founded only on data from the 
optical to the hard X-ray range and, in principle, it could fail 
in describing some specific pulsar properties and the spectral 
behaviour outside this energy range. A useful test is to verify 
how the model reproduces the characteristics of PI described 
in Sect. 5. Fig. 4 shows the model profiles computed at the 
three energies of 2.0, 22 and 85 keV, approximately the mean 
values of the considered energy intervals. The pulse broadening 
is well described and the only deviation is between the phases 
0.025 and 0.075, where the data show a small excess, which 
is s likely due to the simple analytical shape adopted for Cx- 
The search for a better agreement would likely require a more 
complex formula with a greater number of parameters without 
a clear improvement of the general description of the pulsed 
emission. 

A simple approximate way to evaluate the mean photon in- 
dex from the model is to compute the ratio of the summed in- 
tensities of Co and Cx at the two extremes of the energy range: 



/ 



Log(F(E 2 )/F(E l )) 
Log(E 2 /E l ) 



(14) 



Since we used normalised pulse profiles, this spectral index 
would be zero at the centre of PI and therefore we added the 
constants F =1.83 and 2.14 for the MECS and PDS data, re- 
spectively, to match the measured values. In this case the agree- 
ment between the observed points and the model, shown in Fig. 
5, is generally satisfactory, in particular the asymmetric varia- 
tion of F(/) with the softest value at the PI centre. Both these 
results are due to the Cx contribution in the phase interval of 
PI necessary to represent the pedestal apparent in the hard-X 
to low-energy y rays. These results show that observations are 
well reproduced if the spectrum of this pedestal is similar to 
that of Ip, and therefore it is reasonable to consider them all as 
belonging to the same component. 

6.4. The total pulse profile in the MeV range 

Another important verification of the model is the compar- 
ison of pulse profiles and spectra at energies higher than 
the BeppoSAX-NFI range. In particular, we considered the 
low-energy y-ray data obtained by IBIS-PICsIT on board 
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INTEGRAL (Kuiper et al. 2003) and by COMPTEL-CGRO 
(Kuiperetal. 2001). 

We computed the expected pulse shapes from the two com- 
ponent model at the three energies 0.306, 0.830 and 1.90 MeV, 
which are inside the energy ranges of the considered pulse pro- 
files (0.260-0.364), (0.75-1.0) and (1.0-3.0) MeV. We used 
the spectral parameters of Table 4 and the value of p derived 
from Eq. (7). The resulting phase histograms are shown in the 
three panels of Fig. 10 together with the data. We see that the 
model profiles match well the most important features such as 
the heights of the two main peaks and gives an acceptable rep- 
resentation of the Ip region. Differences are generally within 
the observational uncertainties, because of the statistics avail- 
able in these bands. We see therefore that our model is able to 
extrapolate well the pulsed emission up to the MeV range. 

6.5. The emission at optical-UV frequencies 

For a fully consistent description of the Crab pulsed emission 
we evaluated how the spectral laws derived in the previous 
Sections extrapolate in the optical-UV range. At these frequen- 
cies only the contribution expected from Co is important, that 
of Cx being smaller than one order of magnitude. A simple ex- 
trapolation of the model gives a flux lower than that derived 
by the optical-UV observations by Percival et al. (1993) and 
Sollerman et al. (2000). As explained in Sect. 3.2, however, 
reddening effects are very important and even a slight change 
of the extinction parameter can modify the resulting spectrum. 

In order to match the optical points it is not possible to 
reduce the value of b because it would be inconsistent with 
the X-ray data. The simplest way is to assume that at frequen- 
cies lower than ~0. 1 keV the log-parabola tends to approach a 
power law. A possible spectral distribution, including optical- 
UV points, is shown in Fig. 1 1 . 

7. Extension of the two-component model at 
higher energies: the need for other components 

At energies higher than ~30 MeV the model developed in the 
previous section based only on the two components Co and 
Cx is not adequate to describe the emission properties of the 
Crab pulsar. It has been known since early observations that 
the y-ray pulse shape is similar to that of Co, although some 
minor differences are present. The results of COMPTEL and 
EGRET observations (Kuiper et al. 2001) provided high en- 
ergy pulse profiles and spectra good enough to develop a more 
complete model to describe the overall properties of Crab. We 
recall that there are some similarities between the properties 
of the pulsed emission in the hard X-rays and in the high en- 
ergy y-rays. Kanbach (1999) showed that the relative flux of 
P2 to that of PI increases with energy. The pulse profile above 
1 GeV compared to that above 0.5 GeV shows a clear excess 
of Ip and P2 with respect to PI (see, e.g., Fig. 12 of Kuiper et 
al. 2001), Moreover, above ~5 GeV, P2 seems to be the only 
dominant feature in the pulse profile, despite the small number 
of detected events (Thompson 2004). 

To explain these properties in a consistent way, from the 
optical to the GeV band, we extend the previous two compo- 
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Fig. 10. Comparison between the Crab pulse profiles in the en- 
ergy ranges 0.260-0.364 (IBIS-PICsIT data from Kuiper et al. 
2003, upper panel), 0.75-1.0 MeV and 1.0-3.0 MeV (central 
and lower panel, COMPTEL-CGRO data from Kuiper et al. 
2001) computed using the two component model. 



nent model assuming that there are two other high-energy spec- 
tral components, hereafter indicated by Co y and Cx y , which are 
strictly related to Co and Cx, respectively, because their pulse 
shapes are similar and the spectral distributions are approxi- 
mated by log-parabolic laws, but shifted in energy. These ad- 
ditional two components imply at least another six adjustable 
parameters, i.e. the peak energies, curvatures and normaliza- 
tions of the two log-parabolae, but the statistical quality of the 
available data does not allow a good estimate of all of them. 
We assumed, therefore, that their curvatures are equal to that 
of the corresponding low energy components (b = 0.16) and 
found peak energies and normalization factors able to describe 
well the observed spectra. 

We first adapted this four component model to the total 
pulsed spectrum and evaluated the best normalization factors 
of Coy and Cx y - First we noticed that to be consistent with 
the upper limits to the pulsed emission in the TeV range (e.g. 
Lessard et al. 2000, Aharonian et al. 2004) it was necessary 
either to increase the values of b or to include an exponential 
cut-off in the Cx y spectrum. We preferred to follow the latter 
approach, because there is no simple way to fix the curvatures, 
and found that a cut-off energy E c — 15 GeV gives a reasonable 
fit, although this value cannot be precisely evaluated. The same 
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Table 5. Main parameters of the four component model of the 
pulsed emission of the Crab pulsar. 



10 10" 10 10 10 
Energy (keV) 



10" 10' 10° 10" 10" 



Fig. 11. The broad-band spectral energy distribution of the total 
averaged pulse of the Crab Pulsar, with the four components of 
the model. Data points are from BeppoSAX, Compton-CGRO 
and FIGARO II. We also included Te V upper limits (downward 
open triangles) from Aharonian et al. (2004) and the optical 
data (green circles) from Percival et al. (1993), together with 
the best fit after reddening corrections (solid cyan line). Also 
shown (brown dash-dotted line) is the optical dereddened data 
from Sollerman et al. (2000). A possible low frequency extrap- 
olation of the two component SED to match these optical data 
is plotted. 




10 10' 10 

Energy (keV) 



Fig. 12. The broad-band spectral energy distribution of PI 
with the four components of the model. Data points are from 
BeppoSAX, INTEGRAL-ISGRI and Compton-CGRO. 



cut-off energy was also assumed for Co y , although it was not 
required by the data, because the spectral steepening of the log- 
parabola reduces the relevance of this component above a few 
GeV. The resulting spectral energy distribution compared with 
BeppoSAX, COMPTEL and EGRET data is shown in Fig.lTTl 
We also added some data points obtained from the FIGARO II 
experiment (Agrinier et al. 1990, Massaro et al. 1998) to fill the 
gap between 300 keV and 1 MeV. 



Component 


E P 


E c 


LboiC) 




MeV 


GeV 


erg/s 10 36 


Co 


12.2- lfr 3 




2.8 


C x 


178-10" 3 




3.8 


Coy 


300 


<15 


0.97 


Cxy 


2000 


-15 


0.92 



(*) assuming isotropic emission {oj a = An). 



The spectra of PI, Ip, and P2 were computed using the 
phase intervals of Kuiper et al. (2001) and the intensity ratios 
derived from the pulse shapes of Co and Cx- Figs. 1121 113l and 
E]show the resulting spectral energy distribution together with 
the data derived from the literature and our new analysis of 
BeppoSAX and ISGRI-INTEGRAL data. The agreement with 
the data in the PI and P2 ranges is satisfactory, confirming that 
in these intervals the pulse shapes of the Co y and Cx r compo- 
nents are similar to those of Co and Cx, respectively, consis- 
tent with our assumption. The poor statistics of data in the Ip 
region does not allow a good enough verification of the model. 
The SED of Ip has only two components and is poorly known 
above 10 MeV. However, the model expectation for the Ip in- 
tensity was significantly higher than the observational points. 
Consequently, we reduced the normalization of Cx y by a factor 
of about 2 to approximately match the data. This indicates that 
the profile of Cxy in this range must be different from Cx, and 
it does seem to be sharper. When the same change is applied 
to the spectrum of PI we did not find a significant modifica- 
tion, the relative contribution of Cx y being much lower than 
Coy Note in the SEDs of PI and P2 how the mild spectral 
curvature is well reproduced by the superposition of the four 
components and that the positions of the minima are not the 
same: it is around 10 MeV for PI and around 40 MeV for P2. 
At present, however, the data do not allow an accurate estimate 
of this difference. New observations with richer statistics will 
be useful to constrain this feature more precisely, to provide 
further support for our multi-component model and to evaluate 
the component parameters. All the main parameters of the full 
model are given in Table|5] 

This model is not univocally determined by data. It is ap- 
parent from Fig. 12 that in COMPTEL and EGRET data the 
PI spectrum can be fitted by a single power law with the only 
exception being the point at the highest energy where a cut-off 
may be present, as already proposed by Kuiper et al. (2001). In 
the spectral plots evaluated by these authors assuming a power 
law component with a phase independent spectral index, one 
can see that EGRET data for PI are all in excess with respect 
to the best fit, while practically all those of COMPTEL lie un- 
der it. Moreover, large systematic deviations from this best fit 
are apparent in the leading and trailing wings of PI and P2 (see 
also Fierro et al. 1998). We conclude that the model based only 
on components having spectral distributions approximated by a 
log-parabolic law gives a simpler picture of this complex spec- 
tral behaviour. 

Another test of the model is to verify how it reproduces the 
flux ratios P2/P1 and Ip/Pl, frequently considered in the liter- 
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Fig. 13. The broad-band spectral energy distribution of Ip with 
the components of the model. Data points are from BeppoSAX, 
INTEGRAL-ISGRI and Compton-CGRO. 



Fig. 14. The broad-band spectral energy distribution of P2 
with the four components of the model. Data points are from 
BeppoSAX, INTEGRAL-ISGRI and Compton-CGRO. 



ature to describe the pulsed emission of Crab. We computed 
these ratios for the phase intervals adopted by Kuiper et al. 
(2001), and the results are plotted in Figs.[l5]and[l6]for P2/P1 
and Ip/Pl, respectively. The general behaviours of both ratios 
are well reproduced, in particular the sharp decreases between 
1 and 10 MeV and the following slow increase to the GeV 
range. We stress that the sharpness of the decrease depends on 
the curvature and peak energies of log-parabolic spectra and of 
the adopted cut-off. 

We also calculated the bolometric flux of each compo- 
nent by integrating the log-parabolic spectral distribution over 
all energies. This integral can be easily computed analytically 
(Massaro et al. 2004) and the result is: 

AE -(a + bL 0g E) d£ = _ E 2 F(Ep) _ (15) 

Vo 

The luminosity of each of the four components is L — Fboi X 
d 2 aj a , where d is the distance to the Crab and cj a is the solid 
angle occupied by the emission beam. We take d 2 kpc = 
6.2 x 10 21 cm, whereas the value of a> a < An is difficult to 
estimate because it depends on the emission geometry, which 
is not completely known. The resulting values for the isotropic 
case (<±> a - An) are given in Table [5] but they are upper limits 
to the actual values. The total bolometric luminosity, obtained 
by adding the contributions of the four components, is L — 
8.5 x 10 36 erg s _1 , which drops to ~ 6.8 x 10 35 erg s _1 for 
oj a — 1 str. 

8. Discussion 

The development of a detailed physical model for the spectral 
and phase distributions of the broad-band emission from the 
Crab pulsar is a difficult problem. It requires a precise geomet- 
rical definition of the regions inside the magnetosphere where 
the observed radiation originates, and the knowledge of param- 
eters like the orientation angles between the magnetic axis and 



the line of sight to the spin direction. Several models have ap- 
peared in the literature based on either polar cap or outer gap 
geometries. Usually, these models are focused on reproducing 
either the total spectrum or the phase profile and generally they 
are not fully satisfactory in explaining the complex observa- 
tional picture. Moreover, the possibility that the observed fea- 
tures of the pulsed signal can arise from the superposition of 
two or more distinct components is not taken into account. 

We followed another approach and searched for a possi- 
ble interpretation of the Crab signal based on the superposition 
of two or more components that provides a consistent descrip- 
tion of the spectral and phase distributions. MCLM showed that 
a possible explanation of the energy dependence of the pulse 
shape of Crab in the soft to hard X-rays is that we are observing 
two emission components with different phase and energy dis- 
tributions. In that paper we introduced an empiric model, based 
on a collection of BeppoSAX observations, covering the energy 
range from 0. 1 to about 300 ke V, from which it was possible to 
estimate some properties of the main components. Moreover, 
we showed that the X-ray spectrum of PI presents a mild cur- 
vature well fitted by a log-parabola. At energies higher than 
30 MeV, however, this two component model fails to represent 
both pulse profiles and spectra, as observed by EGRET-CGRO 
(Fierro et al. 1998). To take into account y-ray data, Kuiper 
et al. (2001) proposed the existence of three components, two 
of them having log-parabolic spectra, while the third one with 
a power law spectrum and a phase modulated intensity which 
reaches the highest level in correspondence of the two peaks 
and is almost absent in the Ip region. According to our point of 
view, however, this hypothesis is not fully consistent with the 
data, because the flux increases of Ip and P2 are very similar, 
suggesting that their X-ray emission is dominated by the same 
physical mechanism. Moreover, the parameters describing the 
curvatures of their two log-parabolic distributions are very dif- 
ferent, while the analysis of X-ray spectra reported in Sect. 6.2 
indicates that the curvature is rather stable with phase at an 
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intermediate value between them. Finally, the power-law com- 
ponent does not match well the data in some phase intervals at 
y-ray energies where its contribution is dominant. 

To achieve a more consistent scenario, in the present paper 
we propose a model of the pulsed emission from Crab based 
on four components. It is properly a two double-component 
model because each component pair has similar phase distri- 
butions and spectra shifted in energy. At present, our model 
gives an empiric description of the broad-band properties and 
its validity is based on the limited number of assumptions and 
on the resulting capability to obtain a consistent description of 
the observations outside the energy ranges used to evaluate the 
parameters. For instance, the model extrapolates well the pulse 
profiles in the MeV range (see Sect. 6. 4) and the change with 
energy of the P2/P1 and Ip/Pl ratios. We stress that Crab is not 
the only pulsar that shows a behaviour that can be interpreted 
by a multicomponent model. Harding et al. (2002) proposed 
that the X-ray pulsed emission from Vela originates from two 
non-thermal components, one coincindent in phase with the y- 
ray pulse profile and the other one with the optical. 

The physical processes at the origin of these components, 
and the location in the magnetosphere where they occur, must 
be further investigated. The development of a detailed physical 
model of the high energy emission in the Crab magnetosphere 
is beyond the aim of the present paper, however some general 
indications on it can be derived from our conclusions. The same 
phase distribution assumed for each pair of components, such 
as Co and Coy, suggests that their angular pattern and emis- 
sion sites must be coincident or very close, otherwise the aber- 
ration effects would modify the pulse shapes. However, this is 
not necessarily true because on the trailing last open field line, 
aberration and propagation time effects would cancel to form a 
caustic, as shown by Dyks & Rudak (2003). 

According to a widely accepted scenario, primary electrons 
accelerated in a magnetospheric gap emit high energy pho- 
tons which produce e ± pairs against the magnetic field. These 
secondary particles emit synchrotron radiation in the optical 
to MeV energy range. A first possibility to be considered is 
that Co and Cx are synchrotron radiation from two different 
places whereas the corresponding high energy components are 
inverse Compton upscattered photons. A very early inverse 
Compton model for Crab was proposed by Zheleznyakov & 
Shaposhnikov (1972), when the y-ray emission properties were 
known very poorly, and another model was applied to the Vela 
pulsar by Morini (1983). A more recent development is that of 
Cheng and Wei (1995). Their model is based on the outer gap 
geometry and assumes that optical to hard X-ray photons are 
emitted by secondary e ± pairs created outside the accelerating 
region by high-energy primary curvature photons, y-rays are 
then emitted via a synchrotron self-Compton (SSC) process. 
In this case we expect that Compton scattering of hard X-ray 
photons occurs mainly in the Klein-Nishina regime, and there- 
fore the observed high energy cut-off gives an estimate of the 
maximum energy of electrons. The fact that it is in the GeV 
range could be consistent with their origin from magnetic pair 
production. 

There is, however, another possibility for the origin of the 
high energy components. In a polar cap scenario these photons 
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Fig. 15. The ratio between the fluxes of P2 and PI phase re- 
gions (PI: -0.06-0.04; P2: 0.32-0.43), compared to the predic- 
tions of the model. The data points come from various experi- 
ments (Kuiper et al. 2001). The various extrapolations above 1 
GeV correspond to different values of the cut-off energy of the 
C 0y spectrum: 15 GeV (red), 13 GeV (green), 11 GeV (blue) 
and 9 GeV (violet) - from bottom to top. 
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Fig. 16. The ratio between the fluxes of Ip and PI phase regions 
(PI: -0.06-0.04; Ip: 0.14-0.25), compared to the predictions 
of the model. The data points come from various experiments 
(Kuiper et al. 2001). 

could be those emitted by primary electrons via curvature radi- 
ation and not absorbed by magnetic pair production throughout 
the magnetosphere. The attenuation length for this process can 
be approximated by (Erber 1966): 



x(B ± ,E y ) = 2.3 x 10~ 8 — exp 



4B cr 2mc 2 



3B± Ey 



cm, 



(16) 



where B ± is the transverse magnetic field seen by a photon of 
energy E y , and B tT =4.414xl0 13 is the quantum critical field. 
It is easy to verify that a mean free path of the order of 10 4 
cm for 5 GeV photons is obtained in a transverse field B ± =s 
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7.5x 10 8 G and that it increases very rapidly even for decreasing 
E y . This cut-off in the curvature y-ray spectra at GeV energies 
has been verified since early numerical calculations, including 
an accurate evaluation of the absorption coefficient, for a polar 
cap acceleration (Salvati & Massaro 1978, Massaro & Salvati 
1979). 

Although our multi-component model gives a consistent 
picture of the pulsed emission from Crab it is not yet com- 
pletely determined by observational data. An important test 
will be the study of the spectra and pulse profiles at energies 
higher than a few GeV: in particular, we expect that above about 
5 GeV, P2 would be the dominant feature, while a good mea- 
sure of the flux in the Ip region will be very useful to draw the 
phase structure of CV y . High quality data in this energy range 
will be obtained with the LAT telescope on board the GLAST 
mission to be operative next year. Another interesting test of 
the model could be obtained from phase resolved polarization 
measures in the X-ray range. We know that the optical linear 
polarisation of the Ip differs from those of PI and P2 both in 
strength and direction (Smith et al. 1988; Kanbach et al. 2003). 
According to our model the polarization of the P2 X-ray emis- 
sion must be more similar to that of Ip, because of the higher 
contribution of Cx- New generation high sensitivity polarime- 
try for X-ray astronomy, such as that proposed by Costa et al. 
(2001), could be very useful. Finally, the observation of other 
young spin powered pulsars, if their high energy emission is 
similar to Crab, could also be very useful because we expect 
to observe them from different directions and therefore to see 
other pulse shapes, depending on the various combinations of 
the two components. 
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